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ABSTRACT 


Analytic  approximations  to  the  extinction  efficiency,  Qext-,  for  oriented  and  ran¬ 
domly  oriented  circular  infinite  cylinders  based  on  anomalous  diffraction  are  given.  The 
results  are  compared  with  the  exact  code.  These  results  axe  further  generalized  to  ran¬ 
domly  oriented  elliptical  cylinders.  Using  the  formulae,  Qext  can  be  evaluated  over  104 
times  faster  than  with  the  exact  code.  This  approximation  is  valid  for  complex  refractive 
indices  m  =  n  —  ik,  where  1  <  n  <  oo  and  0  <  k  <  3,  aspect  ratios  from  1  to  4  and  modest 
to  large  particle  sizes.  The  accuracy  and  limitations  of  these  formulae  are  discussed. 


RESUME 


Nous  presentons  une  approximation  analytique  de  l’efficacite  d’extinction  Qext  pour 
des  cylindres  infinis  a  orientation  fixe  et  aleatoire.  Cette  approximation  est  comparee  au  cal- 
cul  exact.  Ces  resultats  sont  etendus  aux  cylindres  elliptiques  orientes  aleatoirement.  Grace 
a  cette  formule,  Qext  peut  etre  evalue  104  fois  plus  vite  que  par  les  methodes  anterieures. 
Nous  avons  verifie  cette  formule  pour  des  indices  de  refraction  complexes  m  —  n  —  ik,  ou 
l<n<ooetO<A:<3,  pour  des  aspects  de  1  a  4  et  pour  des  particules  de  taille  moyenne 
a  grosse.  Ces  formules  sont  basees  sur  la  theorie  de  la  diffraction  anomalique.  L’exactitude 
de  ces  formules  ainsi  que  leurs  limites  sont  discutees. 
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EXECUTIVE  SUMMARY 


The  present  and  future  electro-optic  capabilities  of  fire  control  systems  require  or 
will  require  full  spectrum  obscurants  as  an  effective  countermeasure.  It  is  well  known  that 
spherical  particles,  such  as  those  created  in  phosphorous  based  obscurants,  can  obscure 
efficiently  only  in  the  visible  region  of  the  spectrum.  In  the  near  and  far  infrared  as  well  as 
the  microwave  region  of  the  spectrum,  the  use  of  high  aspect,  highly  conductive  cylindrical 
particles  is  mandatory.  For  example,  carbon  fibers  or  iron  whiskers  are  being  considered. 
Furthermore,  reduced  visibilities  due  to  elongated  ice  crystals  in  cirrus  clouds  or  cylindrically 
shaped  hydrosols  will  be  modelled  in  future  realistic  military  simulations. 

The  objective  of  this  work  is  to  reduce  significantly  the  computational  burden  in 
calculating  the  extinction  from  randomly  oriented  cylindrical  (ROC)  particles  that  are  very 
long  in  comparison  with  the  considered  wavelength.  This  allows  for  the  exploration  and 
rapid  computation  of  the  effects  of  varying  size  distributions  and  material  properties  of  the 
ROC  particles  on  the  performance  of  obscurants  and  on  electromagnetic  propagation  in 
general. 


To  date,  the  estimation  of  reduced  visibilities  due  to  aerosols  or  hydrosols  in  military 
simulation  has  been  extremely  crude  or  non-existent.  This  is  in  part  due  to  the  complex  and 
computationally  inefficient  codes  required  to  estimate  the  extinction.  This  work  overcomes 
both  of  these  restrictions.  For  instance,  a  factor  greater  than  ten  thousand  in  computational 
speed  has  been  achieved  with  minimal  loss  in  accuracy.  Furthermore,  fax  less  technical 
knowledge  of  numerical  computation  is  required. 

The  long-term  goal  of  this  work  is  to  alleviate  the  remaining  constraints  in  the  the¬ 
oretical  consideration  of  non-spherical  aerosols  and  obscurants  of  both  natural  and  artificial 
origin.  This  will  not  only  aid  in  finding  better  obscurants  and  allowing  more  accurate  and 
efficient  simulations,  but  also  in  the  possible  identification  and  remote  classification  of  such 
aerosols. 
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1.0  INTRODUCTION 

The  immediate  objective  of  this  work  is  to  reduce  significantly  the  computational 
burden  in  calculating  the  extinction  from  cylindrical  particles.  This  allows  for  the  ex¬ 
ploration  of  the  effects  of  cylindrical  particles  on  the  performance  of  electro-optical  and 
millimeter  wave  systems  in  obscurants,  ice  clouds  and  hydrosols. 

The  current  calculation  with  randomly  oriented  infinite  cylinders  can  require  large 
computation  times,  especially  for  polydispersions,  larger  particle  sizes  and  refractive  indices. 

We  have  previously  presented  (Ref.  1)  an  analytic  approximation  to  Qext  for  ran¬ 
domly  oriented  spheroids.  This  work  applies  to  particles  with  n  >  1  and  k  >  0  for  arbitrary 
sizes  and  aspect  ratios.  The  basic  approach  was  to  orthogonalize  as  much  as  possible  the 
scattering  physics  into  well-defined  regimes.  The  large  particle  regime  was  successfully  mod¬ 
elled  by  anomalous  diffraction  and  edge  (Fock)  terms.  Here  we  will  take  the  same  basic 
approach  but  applied  to  randomly  oriented  cylinders  instead. 

The  report  is  organized  as  follows:  Chapter  2.0  develops  the  extinction  formula. 
Chapter  3.0  contains  comparisons  of  this  approximation  to  the  exact  code  from  IPHASE 
(Ref.  2),  Chapter  4.0  extends  the  anomalous  diffraction  approach  and  explores  the  ram¬ 
ifications,  and  Chapter  5.0  develops  and  discusses  approximations  to  randomly  oriented 
elliptical  cylinder  extinction.  The  final  chapter,  Chapter  6.0,  gives  the  conclusions  and 


limitations. 
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Tliis  work  was  performed  at  DREV  between  October  1993  and  November  1994 
under  PSC  32A,  EO/IR  Protection  of  Land  Vehicles. 

2.0  DEVELOPMENT  OF  THE  EXTINCTION  FORMULA 


The  physics  and  the  structure  of  the  extinction  formula  for  randomly  oriented  cylin¬ 
ders  will  be  developed  in  this  chapter.  Only  medium  to  large  size  parameters  will  be  con¬ 
sidered.  Hence  these  formulae  will  not  apply  in  the  small  particle  or  Rayleigh  region. 

2.1  Anomalous  Diffraction  Term 

The  anomalous  diffraction  approach  was  first  discussed  by  Van  de  Hulst  (Ref.  3) 
for  spheres.  Furthermore,  it  was  extensively  developed  for  oriented  infinite  cylinders  by 
Stephens  (Ref.  4). 

The  anomalous  diffraction  formula  is  derived  by  assuming  that  the  incident  plane 
wave  is  not  significantly  skewed  in  passing  through  the  scattering  object  and  that,  to  first 
order,  the  effect  of  the  scatterer  is  to  delay  locally  the  phase  of  the  wave  and  attenuate 
its  amplitude  (Ref.  3).  The  strict  limit  of  validity  of  the  formula  is  therefore  the  region 
where  (n  —  1)  <<  1.  The  cylinder  is  in  effect  treated  as  a  slit  normal  to  the  incident  wave 
possessing  a  spatially  dependent  phase  and  amplitude.  The  Fraunhofer  pattern  at  infinity 
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is  then  derived  and  Qext  evaluated  from  the  standard  relations.  For  a  cylinder  with  its  axis 
oriented  at  an  angle  6  with  respect  to  the  direction  of  the  incident  radiation,  this  procedure 
leads  immediately  to  the  following  formula  (Ref.  4)  : 


Qad  =  ?rRe[Hi  (a;)  +  t'Ji(w)] 


[1] 


where  u  is  given  by 


a>  =  2 (m  —  l)x/sin  6 , 


[2] 


the  complex  refractive  index  is  m  =  n  —  ik  and  x  is  the  size  parameter.  In  the  above  Hi 
and  Ji  are  the  first  order  Struve  and  Bessel  functions,  respectively.  These  are  well  defined 
in  Ref.  5.  Clearly,  when  the  index  is  real  only  the  first  term  in  [1]  appears. 


For  random  orientation,  the  angle  averaging  gives 


Qad 


£/2  Qad  sin2  0  dO 
J^2  sin2  0  d.0 


=  4Re 


sin2  6  dO  + 


sin2  6  dd 


[3] 

[4] 


The  first  term  in  [4]  can  be  analytically  integrated  as  follows: 

4  f^2  Hi  (p/  sin  0)  sin2  6  dO  =  2  f°° 

Jo  J 1  Zl\/Z  —  1 


[5) 


by  using  the  simple  transform  z  =  1/  sin2  9.  This  integral  appears  as  a  special  case  of 
equation  2.7.4.6  in  Ref.  8.  We  will  use  a  more  general  approach  initially  for  later  purposes. 
We  will  solve  it  by  using  a  general  hypergeometric  form  of  the  above  Weyl  type  integral. 
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From  Ref.  6,  p.  417,  formula  20.5  (2), 


where  G™n  is  the  Meijer  G-function.  The  Struve  function  is  represented  by  the  following 


G-function  (Ref.  7): 


Using  [7]  and  [6]  in  [5]  we  get 


\U-i)d* 


-  2TbaHTUt-h)-  181 

Following  Ref.  7,  p.  230  6.5  (1)  the  above  G-function  can  be  reduced,  giving 

[9] 

This  equation  can  be  directly  obtained  by  equation  2. 7.4.6  of  Ref.  8  as  mentioned  above. 
The  1F2  in  [9]  can  be  further  simplified  by  expressing  it  in  terms  of  Bessel  functions  of  the 
first  kind.  Formulae  to  accomplish  this  can  be  found  (Ref.  8),  giving  finally 

2/~  ^  {*  ”  P~T  [,"(p/2) "  |JoW2)Ji(p/2)  +  (l  -  £)  J?0/2)] } . 


The  above  equation  is  the  full  solution  to  Qad  for  real  indices.  It  is  remarkably  simple  and 
is  readily  computable. 


The  second  term  in  [4]  is  treated  as  follows: 


/•«•/ 2  re 

4i  J  Ji(p/sin^)sin2^d0  =  2  i  J 


Ji  (p^/z)  dz 
z^Jz  —  1 


[11] 
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[12] 


This  integral  in  not  explicitly  in  the  integral  tables,  so  we  will  use  [6]  and  the  G-function 
representation  of  J\: 

Jl(*)  =  ^2  ^"4"|  1  _  -  [13] 

This  gives 


‘fW- f  - 

2  A 

3  11  I  ■ 
2’ 2’  2  / 


1  _1 
2’  2 


dz 


[14] 


This  G-function  can  be  shown  to  be  of  the  logarithmic  type,  which  means  that  no  finite  ex¬ 
pression  in  terms  of  hypergeometric  functions  can  be  obtained.  However,  it  can  be  formally 
expressed  as  an  infinite  series  giving  for  [14]  the  following: 

£  KF  rati  [*** + ^  + + rlr + 2T1] }  m 

In  the  above,  tp(z)  is  the  digamma  function  and  (a)*  is  the  Pochhammer  symbol  (Ref.  5). 
Due  to  the  slow  convergence  of  the  sum  in  [15],  for  even  modestly  large  p,  it  is  not  of  great 
computational  value.  Therefore,  another  simpler  approximation  is  required. 


We  go  back  to  the  physical  level.  When  there  is  no  absorption,  there  is  constructive 
and  destructive  interference.  This  produces  the  familiar  oscillations  in  Qext-  This  is  well 
described  by  [10].  Absorption  affects  the  interference  spatially  by  reducing  the  amplitude 
of  the  partial  waves  inside  the  particle.  This  has  the  effect  of  damping  the  magnitude  of 
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the  oscillations  in  Qext-  For  very  large  absorptions  tbe  particle  is  opaque,  no  waves  are 
transmitted  and  hence,  by  Babinet’s  principle,  Qext  —  2.  To  first  order  we  will  neglect  the 
spatial  dependence  that  absorption  has  on  the  interference.  We  will  do  this  by  averaging 
the  absorption  length  2 ky/ sin  6  spatially  and  over  angles.  Here  y  is  a  cord  length  through 
the  cylinder.  The  angular  average  is  given  by 

So'2  (gSQ  si n20dd  =  8 ky 
Jq/2  sin2  6  d0  * 

Then  doing  the  spatial  averaging  to  get  the  final  averaged  absorption  length  A: 

_  2ffS kfdz_ 

2  giz 


‘ x 2  —  z2  dz 


=  2  kx. 


This  simple  average  gives  e  *  as  the  approximate  damping  factor.  Note  that  averaging 
e-2ky/sm6  pro(juces  another  logarithmic  G-function  of  similar  complexity  as  [15]. 


The  oscillating  terms  can  be  isolated  from  Qac i  (for  real  m)  by  subtracting  from  it 
the  constant  ‘Babinet’  term,  which  is  2.  These  can  then  be  damped  by  e~x  and  then  the 
‘Babinet’  term  added  back.  Thus  the  approximation  to  Qad  including  absorption  becomes 


Qad  »  2  +  e 


{i  -  x  [ JoO>/2)  -  |j0(p/2)  Up! 2)  +  (i  -  j;(p/2)]  }  -  2) 


where  p  =  2(n  —  l)x.  Note  that  p  is  now  always  real. 
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2.2  Edge  Effects 

For  a  particle  whose  typical  size  is  much  larger  than  the  wavelength,  the  edge 
cannot  be  treated  as  sharp  and  the  effect  of  the  curvature  of  the  object  must  be  included. 
Jones  (Ref.  9)  has  shown  how  to  estimate  these  edge  effects  for  convex  bodies.  As  he  clearly 
explains,  near  a  glancing  point,  the  body  will  behave  approximately  like  a  cylinder  with 
its  axis  perpendicular  to  the  surface  normal  (at  the  glancing  point)  and  the  direction  of 
propagation  of  the  incident  wave.  The  scattering  and  absorption  from  all  such  cylindrical 
sections  can  be  integrated  around  the  projection  of  the  penumbral  curve  to  give  the  total 
contribution  of  the  edge  effects  to  the  extinction  cross  section.  The  projection  is  onto  the 
plane  normal  to  the  incident  ray  direction.  The  energy  scattered  per  unit  length  of  the 
cylinder  by  region  on  its  surface  around  the  glancing  point  is  proportional  to  cxeR1^3  for 
the  TE  mode  and  cxmR 1^3  for  the  TM  mode,  where  R  is  the  radius  of  curvature  of  the 
cylinder.  The  constants  cxe  and  ctm  are  the  first  order  Fock  terms.  It  can  be  shown 
that,  for  convex  bodies,  randomly  oriented  or  illuminated  by  a  randomly  polarized  beam, 
(cxe  +  ctm)/2  =  ce,  where  ce  is  a  universal  function  of  refractive  index.  Following  Jones, 
the  edge  contribution  to  the  extinction  efficiency  of  the  body  is 

Q'dae=J  JpR1/3ds  [19] 

where  S  is  the  projected  area  of  the  body  on  a  plane  normal  to  the  direction  of  propagation 
of  the  incident  wave  and  s  is  the  arc  length  along  the  projection  of  the  shadow  boundary 
P  on  that  same  plane. 
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From  the  basic  geometry  of  the  cylinder,  the  terms  in  [19]  become 


S  =  2x1  sin  0,  ds  =  2  sin  ddl,  R 


sin2  0 


where  /  is  a  unit  length  along  the  cylinder  axis.  Thus  Qedge  becomes 

Qedge  —  :  a\2lZ  ’  [21] 

(xsm  Oy/* 

The  angle  averaging  of  [21]  to  obtain  Qcdgc  is  straightforward: 

^  _  f£/2QedgeS  sind  d0  _  2  r(7/6)  ce  1.159595ce 

S  sin  0  dO  ^(5/3)  *2/3  ~  *2/3  '  J 

This  universal  function  ce  is  not  known.  Since  ce  is  independent  of  shape,  this 

could  be  accomplished  by  studying  the  sphere  alone.  This  however,  probably  produces  a 

function  that  is  much  too  complex  to  be  included  in  a  simple  formula.  We  therefore  use  the 
value  for  optically  soft  particles  c*0-^  =  co  =  0.996130,  (Refs.  11-12).  and  (Ref.  12).  An 
added  complication  to  Qedge  arises,  however,  when  we  consider  cylinders  with  small  complex 
phases  xbx.  This  means  that,  for  randomly  oriented  cylinders,  \ib\x  —  2j m  —  l|x  <<  1.  This 
occurs  since  what  we  have  been  calling  an  edge  effect  is  in  fact  the  field  distortion  around 
the  boundaries  of  the  particle,  and  hence  its  behaviour  for  small  phase  is  different  from  that 
for  large  phases.  A  simple  semi-empirical  model  of  the  above  two  behaviours  is: 


^ edge 


1.159595cp 
x2/3  +  x„it 
Pmax 


Here  pmax  is  the  value  of  the  phase  at  the  first  maximum  of  Qext  for  real  indices.  Thus 
xcrit  is  the  size  parameter  that  is  approximately  midway  between  0  and  the  first  peak. 
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Hence  [23]  is  designed  to  have  the  proper  asymptotic  behaviour  for  a;2/3  >>  a ;cr,-t  and  is 
just  a  constant  for  a;2/3  <<  xcr,t.  For  spheres,  pmax  «  4  and  for  randomly  oriented  infinite 
cylinders  pmax  »  Thus  [23]  for  randomly  oriented  infinite  cylinders  becomes 

1.159595cq 
x2/3  +  7r/(4|m  —  1|)’ 

The  large  particle  limit  for  the  randomly  oriented  infinite  cylinders  becomes 

Qext  2  +  ~Qedge-  [25] 

We  now  wish  to  produce  a  term  T  which,  when  it  multiplies  [18],  gives  the  same  limit  as  [25] 
without  diverging  as  the  size  parameter  goes  to  zero.  We  found  that  (Ref.  10)  an  adequate 
expression  for  our  purposes  is 

T  =  2-  e~^<W2.  [26] 

Other  orderings  of  the  angle  averaging  of  the  anomalous  diffraction  and  edge  term  expres¬ 
sions  and  their  products  could  have  been  done,  but  complex  expressions  such  as  [15]  and 
worse  occur  (see  Ref.  1). 
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2-3  Extinction  Formula 

From  the  formulae  in  the  previous  two  sections  the  final  formula  for  Qext  of  ran¬ 
domly  oriented  infinite  cylinders  of  medium  to  large  size  becomes 

[2+e-“*(§^{l-£f  [j?W2)-|W2)Ji(p/2)+(l-£)  Jf (/r>/2)] } -2) ]  (2-e-W*) 

[27] 


This  result  is  valid  for  n  >  1  and  k  >  0. 


It  is  of  interest  to  derive  the  large  p  limit  of  [18]  which  is  the  backbone  of  [27] 
and  thereby  determining  the  peaks  and  valleys  of  the  extinction  curve.  We  do  this  by 
substituting  Hankel’s  asymptotic  expressions  for  the  Bessel  function  equation  9.2.5  (Ref.  5). 
The  result  to  order  1  /p2  is 


Q  cylinder 
ad 


*  acosP  ,  3  A-46inp\ 
2-4  P  +2 


[28] 


This  expression  is  important  since  it  is  similar  to  the  anomalous  diffraction  expression 
derived  by  van  de  Hulst  (Ref.  3)  for  spheres: 


^sphere 


=  2_4?^  +  4 


—  cosp^ 


[29] 


Notice  that  they  are  nearly  the  same  apart  from  an  approximate  phase  shift  of  7r /2  in  p.  This 
implies  that  the  peaks  and  valleys  of  Qext  are  also  shifted  in  the  randomly  oriented  cylinder 
case  by  approximately  tt/2.  Indeed  the  nth  peak  is  at  about  27r(n  —  1/2)  for  cylinders 
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versus  2ir (n  —  1/4)  for  spheres.  Similarly,  the  valleys  are  at  about  2nn  for  cylinders  versus 
2tc (n  +  1/4)  for  spheres. 


3.0  RESULTS 

The  complete  formula,  as  presented  in  the  previous  chapter,  gives  correct  asymp¬ 
totic  behaviour  for  both  large  and  small  |n  —  1|,  k  <  1  and  x  >  1.  In  studying  the  error 
behaviour  of  the  approximation,  the  mid-ranges  of  j m  —  1|  and  x  are  of  greatest  interest.  In 
this  chapter  we  will  compare  the  analytic  approximation  with  the  exact  randomly  oriented 
cylinder  code  as  implemented  in  IPHASE  (Ref.  2). 

Figure  1  shows  the  comparison  of  Qext  versus  x  for  a  refractive  index  of  n  =  1.3,  close 
to  that  of  ice  in  the  visible.  It  is  clear  that  the  error  decreases  for  large  x.  The  largest  errors 
are  near  and  around  the  first  two  peaks.  More  of  the  scattering  physics  must  be  considered. 
All  the  peaks/ valleys,  while  they  are  well  modelled,  are  underestimated/overestimated.  This 
occurs  since  anomalous  diffraction  assumes  that  the  rays  are  not  refracted  by  the  particle 
and  therefore  overestimates  the  phase  difference  as  a  function  of  angle.  In  the  next  section, 
the  anomalous  diffraction  formulation  will  be  extended  to  account  for  this  effect. 

For  small  k  and  large  n  one  can  expect  a  significant  contribution  from  surface  waves. 
This  effect  is  shown  in  Fig.  2  for  an  index  of  m  =  1.8—0  i.  Unlike  those  for  a  sphere  (Ref.  1), 
these  resonant  waves  do  not  contribute  as  much  and  hence  the  present  approximation  does 


much  better. 
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FIGURE  1  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.3.  Qtxt  of  ice  in  the  visible. 

For  large  n,  body  resonances  can  occur.  These  are  sometimes  called  morphology 
dependent  resonances  (MDR).  For  spheres,  these  occur  near  nx  =  lw,  where  l  is  a  natural 
number.  For  cylinders,  in  any  orientation,  they  occur  at  the  first  roots  of  the  integer  order 
Bessel  functions.  Hence,  the  MDR’s  occur  when  Ji(nx)  «  0.  For  example  the  first  MDR 
occurs  when  nx  «  2.40. 

Figure  3  shows  an  example  of  an  incipient  MDR  on  the  first  diffraction  peak  of 
Qext •  Here,  m  =  3  —  Oi.  Note  that,  despite  the  significant  perturbation  in  the  transition 
region  due  to  the  MDR,  the  approximation  is  excellent. 

The  next  two  figures  are  for  ice  cylinders  in  the  3  —  5 /mi  and  8—12  /z m  regions, 
respectively.  In  Fig.  4  the  diffraction  peaks,  as  in  the  previous  figures,  are  underestimated. 
Figure  5  is  a  case  with  large  absorption  and  shows  that  the  simple  exponential  damping 


Qext 
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FIGURE  2  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.8  and  an  aspect  ratio  of  1;  significant  surface 
waves. 


FIGURE  3  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  3  —  Oi;  incipient  MDR  at  x  «  0.8. 
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FIGURE  4  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.45  —  0.0015  i.  Qext  of  ice  in  the  mid-infrared. 


FIGURE  5  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.10  —  0.134  i.  Qext  of  ice  in  the  far-infrared. 
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term  e~2kx,  in  [18],  is  quite  successful  in  modelling  the  rise  of  the  extinction. 

Figure  6  represents  the  calculation  for  typical  hydrosols.  Many  particles  that  are 
suspended  in  water  have  an  absolute  index  of  about  1.4.  Thus  the  index  to  be  considered 
in  the  scattering  process  is  relative  to  water.  Hence,  in  this  example,  m  fa  1.4/1.33  «  1.05. 
This  is  the  region  where  anomalous  diffraction  applies  and  the  excellent  agreement  with  the 
exact  results  in  Fig.  6  testifies  to  this. 

The  last  two  figures  in  this  section  are  for  carbon  fibers  in  the  ultraviolet  to  visible 
and  3  —  5  /zm  regions.  Figure  7  has  even  higher  absorption  than  that  of  Fig.  5  but  is  still 
well  modelled.  As  in  all  the  above  cases  for  x  <  1  we  do  not  expect  very  good  agreement 
since  there  has  been  no  attempt  to  model  the  Rayleigh  region.  In  Fig.  8,  the  absorption 
is  extremely  high  and  is  high  enough  for  the  particle  to  be  reflective.  The  overestimate  of 
Qext  for  x  <  3  is  due  to  a  breakdown  of  our  model  of  the  edge  effect. 

All  our  approximate  Qext  diagrams  in  this  chapter  were  produced  at  a  rate  of  greater 
than  104  times  faster  than  by  the  exact  code.  The  execution  time  for  the  approximation 
on  an  Intel  i860  40  MHz  coprocessor  is  about  70  fis  and  is  independent  of  particle  size  and 
index.  This  large  speed  up  factor  is  essential  if  a  complex  computer  model  requires  a  large 
number  of  computations  such  as  in  a  realistic  military  simulation. 
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FIGURE  6  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.05.  Q^t  of  hydrosols  in  the  visible. 


FIGURE  7  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  1.57  —  .98  i.  Qext  of  carbon  fibers  in  the  ultra¬ 
violet  to  visible. 
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3.5  |j 

3  M=4.0— 2.9i 


SIZE  PARAMETER 

FIGURE  8  -  Comparison  between  approximation  and  randomly  oriented  cylinder 
code  for  an  index  of  4  —  2.9  i.  Qcxt  of  carbon  fibers  in  the  mid-infrared. 


In  the  previous  chapter  it  was  pointed  out  that  the  standard  anomalous  diffraction 
underestimates/overestimates  the  peaks/valleys  in  the  extinction  curve  since  only  unde¬ 
viated  rays  are  considered.  In  this  chapter  the  extended  anomalous  formalism  will  be 
applied  first  to  oriented  infinite  cylinders  and  then  the  orientation  will  be  randomized.  This 
approach  was  first  used  to  improve  the  extinction  formula  for  randomly  oriented  prolate 
spheroids  (Ref.  10).  This  work  also  underlines  the  fact  that  ignoring  the  deviated  ray 
produces  unacceptable  errors  for  convex  elongated  oriented  particles. 
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4.1  Oriented  Infinite  Cylinders 

The  extended  anomalous  diffraction  formula  for  oriented  infinite  cylinders  is  given 
by  [1],  except  that  the  phase  difference  u>  now  takes  into  account  the  refraction  of  the  central 
ray  through  the  cylinder: 

Qad  =  7rRe[Hi(u>)  +  iJi(w)]  [30] 

where  u>  is  now  given  by 

u>  =  tpx  =  2x[(m?  —  cos2  0)1/2  —  sin  0\.  [31] 


Using  [21]  and  [23],  the  edge  term  for  oriented  infinite  cylinders  in  the  extended 
anomalous  approximation  is 


'Qedge 

Xcrit 


CO 


(a;2/3  _|_  Xcrit)  sin2/3  6 

_ 3I6 _ 

4[(m2  —  cos2  0)1/2  —  sin  6\ ' 


[32] 


Since  the  first  maximum  of  the  Hi (2)  occurs  at  about  z  =  3.6,  the  above  has  pmax  =  3.6. 


To  be  complete,  [31]  for  the  standard  anomalous  approach  is 


10  =  =  2  x(m  —  1)/  sin  6. 


[33] 


and  the  edge  term  for  undeviated  rays  is  [32]  but  now 


Xcrit  = 


3.6  sin  6 
4|m  —  1| 


[34] 
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FIGURE  9  -  Effect  of  extended  anomalous  diffraction  on  Qext-  The  index  is  1.3  and 
orientation  angle  is  90°. 


The  effects  of  this  extension  on  the  extinction  for  a  particular  orientation  can  be 
dramatic,  as  shown  in  the  next  sequence  of  three  figures.  The  orientation  angle  is  decreased 
from  90°  to  30°.  For  Fig.  9  there  is  no  difference  between  the  standard  and  extended 
anomalous  diffraction  approaches  since  there  is  no  deviation  of  the  central  ray.  For  60°, 
Fig.  10,  the  phase  of  the  standard  anomalous  diffraction  solution  is  clearly  incorrect  and 
this  error  accumulates  at  larger  x.  At  30°,  Fig.  11,  the  situation  worsens.  This  continues 
for  smaller  orientation  angles. 
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FIGURE  10  -  Effect  of  extended  anomalous  diffraction  on  Qext-  The  index  is  1.3  and 
orientation  angle  is  60°. 


FIGURE  11  -  Effect  of  extended  anomalous  diffraction  on  Qext  •  The  index  is  1.3  and 
orientation  angle  is  30°. 
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4.2  Random  Orientation 

In  the  angle  averaged  case,  the  effects  demonstrated  by  Figs.  9-11  can  still  be  seen. 
This  is  shown  in  Figs.  12  and  13  for  m  =  1.3  and  m  =  1.8,  respectively.  Taking  account 
of  the  deviated  ray  improves  both  the  shape  and  amplitude  of  the  estimate  Qext •  This  is 
more  evident  at  the  higher  index. 

Since  the  Jacobian,  sin2  9,  weights  towards  90°,  the  effect  of  ray  deviation  is  not 
as  dramatic  as  in  the  oriented  case.  Additionally,  the  integral  produced  by  the  angle  av¬ 
eraging  does  not  have  an  exact  analytic  solution.  An  approximate  analytic  formula  has 
been  obtained  (Refs.  13-14),  but  is  considerably  more  complex  than  [27]  and  has  an  error 
comparable  to  the  improvement.  Hence,  it  is  satisfactory  to  consider  only  undeviated  rays 
in  approximating  Qext  for  random  orientations  if  a  simple  efficient  formula  is  of  paramount 


importance. 
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FIGURE  12  -  Effect  of  extended  anomalous  diffraction  on  Qext ■  The  index  is  1.3, 
random  orientation. 


5.0  n 


FIGURE  13  -  Effect  of  extended  anomalous  diffraction  on  Qexu  The  index  is  1.8, 
random  orientation. 
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5.0  ELLIPTIC  INFINITE  CYLINDERS:  AN  EXTENSION 

The  power  of  the  anomalous  diffraction  approach  can  be  extended  in  a  straightfor¬ 
ward  manner  to  infinite  cylinders  of  elliptical  cross  sections.  Since  the  extended  anomalous 
diffraction  generally  leads  to  complicated  expressions  (Ref.  1),  we  will  present  only  formula 
for  the  standard  anomalous  diffraction.  We  will  first  set  up  the  formalism  for  the  oriented 
case  and  then  do  the  angle  orientation. 


5.1  Oriented  Case 


Let  the  infinite  axis  of  the  elliptic  cylinder  be  oriented  at  an  angle  9  with  respect 
to  the  direction  of  the  incident  radiation,  as  was  the  case  for  the  circular  cylinders.  Take 
a  cross  section  normal  to  the  infinite  axis  of  the  cylinder.  Let  </>  be  the  angle  between 
the  semi-major  axis  of  the  elliptical  cross  section  of  the  cylinder  and  the  plane  defined  by 
the  direction  of  the  incident  radiation  and  the  infinite  axis  of  the  cylinder.  Let  p  be  the 
projection  operator  of  this  elliptical  cross  section  onto  the  shadow  plane  (Ref.  3).  Then, 


P  = 


\j cos 


2  <f>  r2  sin2 


<t> 


[35] 


where  r  is  the  ratio  of  the  semi-major  axis  size  parameter  a  to  the  semi-minor  axis  size 
parameter  b  of  the  elliptic  cross  section.  The  anomalous  diffraction  expression  Qa(i  for  an 
infinite  elliptical  cylinder  oriented  at  angles  6  and  <j>  is 


Qad  —  5rRe[Hi(a)  +  iJi(o:)] 


[36] 
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where  a  is  given  by 


_  2(m  —  1  )rb 
a  p&m  6  1 


[37] 


Thus  Qad  is  the  same  as  for  the  circular  cylinder  but  with  x  replaced  by  rb/p.  The  edge 


term  is  simply  derived  to  be 


These  equations  allow  us  to  study  the  randomly  oriented  case. 


5.2  Random  Orientation 


The  angle  averaging  must  be  done  over  two  angles,  6  and  <f>.  It  is  evident  that  the 
averaging  over  0  will  give  [27]  with  p  — ►  2 (n  —  l)rb/(p sin#),  x  —*•  rb/p  and 

1.159595  cor2/3 


Q  edge 


p2[&2/3  +  T/(4\m  -!|)]' 


[39] 


This  leaves  the  angle  averaging  over  <f>  to  be  performed.  An  analytic  expression  can  be 
found,  but  is  complicated.  However,  it  can  be  numerically  integrated  or  approximated  in 
various  limits.  The  latter  will  be  discussed  in  the  next  section. 


Figures  14  and  15  are  computations  of  the  extinction  efficiency  predicted  by  stand¬ 
ard  anomalous  diffraction  plus  edge  term  for  randomly  oriented  elliptic  cylinders.  Figure  14 
compares  the  results  between  a  circular  cylinder  and  an  elliptic  cylinder  with  r  =  2.  Figure 
15  is  the  same,  but  r  =  3  for  the  elliptic  cylinder.  Note  that  the  peaks  are  shifted  towards 
smaller  b  and  the  amplitudes  are  reduced  with  increasing  r.  Most  of  the  peak  shift  occurs 


between  r  =  1  and  r  =  2. 
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4.0  u 

J 


FIGURE  14  -  Comparison  between  randomly  oriented  circular  cylinders  and  elliptic 
cylinders.  The  index  is  1.3  and  r  =  2. 


4.o  q 
] 


FIGURE  15  -  Comparison  between  randomly  oriented  circular  cylinders  and  elliptic 
cylinders.  The  index  is  1.3  and  r  =  3. 
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6.3  Approximations 


If  we  make  the  same  substitution,  as  was  done  in  the  previous  section  but  into  the 
asymptotic  form  [28]  and  then  perform  an  approximate  angle  averaging  over  the  first  terms 
only  one  obtains  a  simple  and  reasonably  accurate  expression  for  the  extinction  efficiency. 
Thus  the  completely  angle  averaged  extinction  efficiency,  <  Qad  >,  is  approximately 


_  ^4  cos (atr/p)  2  , ,  r^, 

<  Qad  >«  2  -  ^  Jo  —  P  dcj>  [40] 

where  E(e2)  is  the  complete  elliptic  integral  of  the  second  kind,  e2  =  1  —  1/r2  is  the  square 
of  the  eccentricity  and  a  =  2 (n  —  1)6. 


The  integral  in  the  second  term  of  [40]  must  be  approximated.  By  rewritting  p  = 
yj (cos2(<£)  +  r2  sin2(<jS)  as  r-i/1  —  e2cos2((f>)  =  rp'  we  can  expand  1/p  in  terms  of  e2,  which 
is  always  between  0  and  1  and  hence  small.  Thus 


i«i[l +  /(£*)  «**(,«]. 


[41] 


If  the  Taylor  expansion  is  used  then  /(e2)  =  e2/2.  This  function,  as  will  be  shown  later, 
can  be  better  optimized. 

Putting  [41]  into  the  integral  in  [40],  we  get 

f  cos(ar/p)p 2  d(f)  &  —  f  cos[a(l  +  f(e2)  cos2(0))][l  —  e2  cos2(<£)]  d<j>.  [42] 

j o  glt*  a  j q 


After  some  significant  algebra  the  evaluation  of  [42]  requires  performing  the  follow- 
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ing  integrals: 

I\  =  f  cos(2  (3  cos  2(<£))  d4>  =  —  cos(/3)  J0(/3) 

Jo  2 

Jrfr  /2 

’  sin(2/3  cos2(<£))  d<f>  =  -  sin(/3)  J0(/3) 

o  ^ 

/3  =  JQ  1  cos(2/?cos2(<£))  cos2(<£)  #  =  ^2^3  1/2  1^2  ~  ^  ' 
=  ^[cos(/?)  J0(/3)  -  sin(/3)  Ji(/?)] 

Ti  ~  1  sin(2/3cos2(^))  cos2 (<j>)  dcfr  =  ^2-Ps  3/2^3/2^2  “  ^ 

=  -J[sin(/3)  Jo(^)  +  cos(/3)  Jx(/3)] 


[43] 

[44] 

[45] 

[46] 


In  the  above,  Jx  to  J4  can  be  obtained  by  using  Mathematica  V.2.2.  The  reduction 
of  the  two  2-F3  hypergeometric  functions  in  I3  and  I4  requires  breaking  them  into  their  even 
and  odd  parts  (MacRobert’s  theorem)  and  then  using  the  reduction  tables  (Ref.  8). 

Putting  these  results  into  [40],  we  finally  obtain 

<  Qad >  *  (i-|)  cos(/S) +  m 

_si^2)  ^1  _  Sn(J))J0(f 3)  -  t  co.(0J,(j8)  }  [48] 

where  ,5  =  af(e2). 

There  remains  to  determine  an  optimal  form  for  /(e2).  A  simple  but  effective 
approach  is  to  minimize  the  least-square  error,  £,  of  1/p  or  1  /p'  and  [41]  i.e: 

S  =  Min  |  J  1  [(1  +  /(e2)  cos2(<£))(p')2  -  p'] 2  J .  [49] 
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Thus 

ftp  r  i 

■qJ  =  2  jf  [(1  +  /(e2)  cos2(«£))(p')2  ~  p']  cos2(^>)  (p')2  =  0.  [50] 

Solving  for  /  we  get 

f’r/2  «r/2  /■7r/2 

'  (p'f  COS2 (4>)d(f>-  I  (p')4cos2(<j>)d(j)  /  [  (p1)4  cos4(<f>)  d<j>  =  [J5  -  J6]//7 

0  JO  J  Jo 

[51] 

The  above  three  integrals  I5  to  I7  can  be  solved  again  by  using  Mathematica  and  the  Gauss 
contiguity  relations  (Ref.5).  Thus 

h  -  jf  (p'f'COS2(4>)d4,=  ~  2^1  3/^>3/2  ^ 

=  5?K3  -  ^  +  4e<«<2)  -  <3  -  1&2  +  *W)1 
“  2(4^3)  [2  +  (If  "  3)  '2  +  £1  1521 

h  =  f  (p')4cos2(»d<£ 

Jo 

=  §  [2  -  3^  +  |£-]  [53] 

h  =  /  (pO4  cos4(<^)  d4> 

Jo 

In  the  above  K(e2)  is  the  complete  elliptic  integral  of  the  first  kind.  The  approximation  in 
[52]  is  a  rational  approximation  with  matched  asymptotes  to  simplify  the  computation,  if 
required.  This  gives  a  rational  form  for  f(e2)  that  is  very  close  to  the  Taylor  solution  for 
e  — >  0.  It  is  nearly  a  factor  of  2  larger  than  the  Taylor,  however,  when  c  — ►  1  and  hence  it 


is  where  one  expects  the  greatest  improvement  in  ^  Q  ad  ^  • 
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A  second  approximation,  that  is  good  for  small  e2,  can  be  obtained  by  a  Taylor 
expansion  of  <  Qad  >  about  e2  =  0,  i.e. 

jt  d  <  Qa.d{e2)  >  1 


<  Qad(e2)  >=<  Qad( o)  >  +e2 


dc2 


+  •  •  •  • 


[55] 


e2=0 


This  expansion  will  be  not  only  restricted  to  small  e2  but  also  small  b  since  the  phase 
information  in  <  Qa<i  >  is  also  being  approximated.  With  some  algebraic  manipulation  the 
second  term  on  the  right-hand  side  becomes 


9  <  Qadje2)  > 

de 2 


a 2  7 re*3 


£*= 0 


16 
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F2 


1/2 

2,3 


cr 

T 


a  _ 

—  - lFo 

48  1  2 


3/2 

3,4 


<  Qad( 0)  > 

[56] 


Hence,  with  additional  work,  the  approximation  to  <  Qad(e2)  >  for  small  e2  and  b  becomes 

<  Qad{<?)  >  «  (i  +  J^j  <  Qad{ 0)  >  Jl  (|) 

=  H)  ¥  {>-f  [■*  (I)-!*  (I)  *  (fM1-?)  *(f)]} 

[57] 

An  additional  approximation  can  be  obtained  by  taking  a  formal  asymptotic  expansion  of 
[28]  (Ref.  16).  Following  the  formula  (6.3.28)  and  (6.3.38)  of  (Ref.  16)  we  get 

2 


<  Qad  >  ~  2 


a2E( 


aE(e 2 


[(1  - €,)V4  cos  (v!=7  -  ?)  + cos  (“ ' +  ?)]  + 


k (£2)- 


3\/S  [(1  -  ^ ' sin  (^=7  '  I)  + sin  (“ + j)] ) 


[58] 


This  formula  is  asymptotic  in  a  and  e.  It  is  clear  that,  as  e  goes  to  0,  [58]  does  not  go  to 
the  asymptotic  circular  cylinder  formula  [28].  This  can  be  easily  corrected  by  multiplying 
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FIGURE  16  -  Comparison  between  numerically  integrated  <  Qad  >  and  [58].  The 
index  is  1.5  and  r  =  1.4.  (Note  that  a  =  b ). 


the  relevant  terms  in  [58]  by  1  —  e-vVe2a/2.  ^he  result  is 


<  Qad  >' 


2  + 


a2E(e2) 


(2  -  E2)^)  -  T(£2) 


[(1  - ^  - -  7)  ■ +  -  (“ + !)] + 

|  [(1  -  c2)8*  sin  (^=5  -  f)  +  **  («+!)]}  [59] 


This  formula  is  now  asymptotic  in  a,  for  any  e. 


Figure  16  is  a  comparison  between  the  numerically  calculated  <  Qa<i  >  and  [59]. 
The  numerical  integration  was  carried  out  by  using  a  256-point  Gaussian  quadrature.  It  is 
surprising  that  such  a  simple  expression  as  [59]  can  capture  all  of  the  subtle  features  of  the 


oscillations. 
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This  formula  makes  it  obvious  that  the  extinction  peaks  for  large  aspect  cylinders 
occur  near  tt (n  —  1/4).  This  formula  forms  an  excellent  backbone  to  <  Qext  >.  All  that 
remains  is  to  add  the  effects  of  absorption  and  the  edge  terms  in  a  manner  similar  to  that 
done  for  circular  cylinders.  Hence,  the  addition  of  absorption  is  performed  as  in  [18].  All 
that  is  required  is  the  following  simple  integration  to  obtain  the  averaged  absorption  length 
A: 


The  edge  term  also  requires  just  a  simple  integration.  It  is  set  up  by  using  [39]  and 
changing  the  Xcrit,  7r/(4|m  —  1|)  by  the  semi-empirical  form  7r(l  —  e2/4)/(4|m  —  1|6).  See 
the  disscusion  concerning  [23]. 


<‘-  Qedge 


1.159595  co  f*/2  f  1  \  , 

r4/3[fe2/3  +  T(!  _  €2/4)/(4|m  _  l|)]£(e2)  J0  ^(p/)2  )  P  d<? 
1.159595  co  K(e2) 

r4/3[62/3  +  ff(l  _  e2/4)/(4jm  _  !|)]  jB(e2) 


[61] 


Thus,  the  asymptotic  formula  for  <  Qext  >  is 

<  Qext  >-  [2  -  e-A(<  Qad  >  -2)]  (1  -  e-<G«ds*>/2)  [62] 


Unfortunately,  the  accuracy  of  the  above  formula  cannot  at  present  be  determined 
since  there  are  no  known  results  for  randomly  oriented  elliptic  cylinders.  The  results  to 
date  are  for  a  particular  orientation  only  (see  Ref.  17  and  references  therein). 
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6.0  CONCLUSIONS  AND  LIMITATIONS 

We  have  presented  an  analytic  approximation  to  Qtxt  for  oriented  and  randomly 
oriented  infinite  cylinders  of  circular  cross  section.  In  addition,  a  formula  requiring  only 
one  numerical  integration  is  given  for  randomly  oriented  infinite  cylinders  of  elliptic  cross 
section.  Various  asymptotic  expressions  for  these  formulae  are  also  developed. 

The  large  size  regime  is  modelled  by  anomalous  diffraction  plus  ‘edge’  or  Fock 
terms  as  in  Ref.  1.  The  small  particle  size  limit,  although  roughly  modelled,  is  not  treated 
extensively  since  the  resulting  formulae  would  be  too  complex.  Furthermore,  at  least  in  the 
case  of  circular  infinite  cylinders,  the  exact  code  is  efficient  in  this  regime. 

An  extension  to  the  anomalous  diffraction  approach,  similar  to  that  in  Ref.  1,  is 
also  developed  for  the  case  of  circular  cylinders.  This  demonstrates  that  this  extended 
anomalous  diffraction  is  significant  only  when  the  cylinders  are  oriented. 

The  circular  cylinder  formula  gives  good  results  with  little  loss  in  accuracy  for 
medium  to  large  size  parameters  for  n  >  1  and  0  <  k  <  3  and,  for  elliptic  cylinders,  for 
modest  aspect  ratios  (r  <  4).  Typical  errors  of  about  5%  or  less  are  encountered  in  the 
circular  cylinder  case.  The  errors  observed  have  the  same  behaviour  as  previous  seen  for 
randomly  oriented  spheroids  (Ref.  1).  Errors  for  the  elliptic  cylinder  cannot  be  determined 
at  this  date  since  there  are  no  known  exact  results  (Ref.17).  However,  there  is  no  physical 
reason  we  know  of  that  would  make  the  errors  in  the  formulae  behave  differently  for  elliptical 
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cylinders  than  that  of  either  the  circular  cylinders  or  spheroids. 

If  high  precision  is  not  required,  the  circular  cylinder  formula  is  far  more  economical 
in  computer  time  than  the  exact  random  infinite  cylinder  code.  This  is  even  more  true 
when  considering  elliptical  cylinders  or  poly  dispersions.  Taking  account  of  both  the  range 
of  demonstrated  validity  and  the  accuracy,  the  circular  cylinder  formula  is  superior  to  all 
other  approximations  known  by  the  authors. 

An  even  simpler  form  for  large  optical  sizes  is  also  derived  for  both  the  circular  and 
elliptic  cylinders  which  are  analogues  of  Van  de  Hulst’s  famous  formula  for  spheres.  This 
allows  for  the  simple  determination  of  the  maxima  and  minima  of  the  extinction. 

The  results  of  this  work  demonstrates,  once  again,  the  effectiveness  of  the  combined 
approaches  of  anomalous  diffraction  and  edge  effect  terms  in  modelling  the  extinction  from 
smooth  convex  bodies.  This  is  a  positive  signal  towards  further  work  on  other  shapes  and 
scattering  properties  to  fully  cover  the  broad  spectrum  of  natural  and  artifical  scatterers. 
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